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Abstract 

Two-photon and two-gluon annihilation rates of P-wave scalar charmonium and bottomonium 
up to third radial excited states are estimated in the relativistic Salpeter method. We solved the 
full Salpeter equation with a well defined relativistic wave function and calculated the transition 
amplitude using the Mandelstam formalism. Our model dependent estimates for the decay widths: 
r(x c o -» 2 7 ) = 3.78 keV, T{ X ' c0 - 2 7 ) = 3.51 keV, F( X bo -» 2 7 ) = 48.8 eV and r( X ' 60 - 2 7 ) = 50.3 
eV. We also give estimates of total widths by the two-gluon decay rates: T to t(Xco) = 10.3 MeV, 
r toi (Xco) = 9-61 MeV, r tot (x6o) = 0.887 MeV and T tot ( X ' b0 ) = 0.914 MeV. 
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I. INTRODUCTION 



ation rate, we have 
[13|. First was the 



It is well known that two-photon or two-gluon annihilation rate of heavy quarkonium 
cc or bb is related to the wave function, so this process will be helpful to understand the 
formalism of inter-quark interactions, and can be a sensitive test of the potential model l| . 
With a replacement of the photons by gluons, the finial state becomes two gluon state, which 
will be helpful to give the information on total width of the corresponding quarkonium. 

In previous letter [2J, two-photon and two-gluon annihilation rates of h pseudoscalar 
cc and bb states are computed in the relativistic Salpeter method, good agreement of our 
predictions with other theoretical calculations and available experimental data is found. We 
also found the relativistic corrections are important and can not be ignored. In this letter, 
we extend our previous analysis to include the P-wave ++ scalar cc and bb states, present 
a relativistic calculation of these states decaying into two photons or two gluons. 

For the theoretical estimates of the two-photon or two-gluon annihi 

various methods readily available in hand {3, 4], |], Q, E, Q, [9, [uj, [ill, fl2 

non-relativistic calculation, the corresponding decay width is related to the derivative of 

the non-relativistic P-wave function at the origin, this method will cause large uncertainty 

because in a full relativistic calculation the decay width is related to the full behavior of 
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P-wave function which can be seen in this letter or in Ref. [2J. So the relativistic corrections 
is very important. In recent years, many authors try to focus on the relativistic corrections 
and there are already some versions of relativistic calculation, and they give improved results 
over the non-relativistic methods. In this letter, we give yet another relativistic calculation 
by the instantaneous Bethe-Salpeter method 14|, which is a full relativistic method 15] 
with a well defined relativistic form of wave function. 

In previous letter j^j], we have pointed out that there are two sources of relativistic correc- 
tions; one is the correction in relativistic kinematics which appears in the decay amplitudes 
through a well defined form of relativistic wave function (related to the full behavior of wave 
function, not merely related to the derivative of the wave function at origin); the other rel- 
ativistic correction comes via the relativistic inter-quark dynamics, which requires not only 
a well defined relativistic wave function but also a good relativistic formalism to describe 
the interactions among quarks. The Bethe-Salpeter equation and its instantaneous version, 
Salpeter equation, are well-known tools to describe relativistic bound states. In this letter 
we will solve the full Salpeter equation for ++ state, and use the full Salpeter wave function 
to estimate the annihilation decay width of quarkonium. 
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The form of the wave function is also important in the calculation, since the corrections 
of the relativistic kinetics come mainly through it. We begin from the quantum field theory, 
analyze the parity and charge conjugation of bound state, and give a formula for the wave 
function that is in a relativistic form with definite parity and charge conjugation symmetry. 
Another important thing is how to use the relativistic wave function of bound state to obtain 
a relativistic transition amplitude, since a non- relativistic transition amplitude even with a 
relativistic wave function will lose the benefit of relativistic effects caused by a relativistic 
wave function. The Mandelstam formalism is well suited for the computation of relativistic 
transition amplitude, and we begin with this formulism to give a formula of the transition 
amplitude. 

In Sec. II, we give the transition amplitude in Mandelstam formalism and corresponding 
wave function with a well defined relativistic form. In Sec. Ill, the full Salpeter equation is 
solved, and the mass spectra and numerical value of wave function are obtained. Then the 
two-photon decay width and full width of heavy ++ quarkonium are estimated. In Sec. Ill, 
short discussions and a summary are also given. 

II. THEORETICAL DETAILS 

According to the Mandelstam [16] formalism, the relativistic transition amplitude of a 
quarkonium decaying into two photons (see figure 1) can be written as: 



T 27 = iV3 (iee q ) 2 ' Jj^J tr | x(o) 



^S(pi - k x ) fa + faS(pi - k 2 ) 



(1) 



(27T) 

where kx, k 2 ; £x, £2 are the momenta and polarization vectors of photons; e q = | for charm 
quark and e q — | for bottom quark; p x and p 2 are the momentum of constitute quark and 
antiquark; xio) is the quarkonium Bethe-Salpeter wave function with the total momentum 
P and relative momentum q, related by 

mi 

p 1 = a x P + q, ax ~ 



p 2 = a 2 P - q, a 2 = 



nil + m 2 
m 2 



nix + ^2 

where mi = m 2 is the constitute quark mass of c or b. 



Since p w 
mass system 



I 'on = M, the approximation p w = p 20 = -y is a good choice for the equal 
Having this approximation, we can perform the integration over q 
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to reduce the expression, with the notation of Salpeter wave function ^f(q) = J $ Q -x(?), to 



T 27 = V3(e eg ) 2 /^tr(^ 



1 1 

+ 



i>\- h - mi rl 1 i>x- % - rni r ' \ 
With this relativistic amplitude, the two photon decay width can be written as 



(2) 



T 2 

The general form for the relativistic wave function of scalar state J PC = ++ can be 
written as 8 terms constructed by momentum P, q and Dirac matrix 7, because of the 
approximation of instantaneous, 4 terms with P ■ q become zero, with further constraint 
from Salpeter equation, the relativistic Salpeter wave function ^(q) for ++ state with a 
definite parity (+) and charge conjugation (+) can be written as: 

V(Q) = <Pi(q)£+<P2(q)lojf-—<Pi(q)- (4) 

The wave function <f\{q), ¥2(0) an d bound state mass M can be obtained by solving the 
full Salpeter equation with the constituent quark mass as input, and they should satisfy the 
normalization condition: 

J {ZTTJ^ Till 



where uo\ = \jm{ + q 2 . 

The two gluon decay width of quarkonium can be easily obtained from the two photon 
decay width, with a simple replacement in the photon decay width formula 

9 

4 2 ^ 2 tr\ 

e q a — ► - a s . (6) 



III. NUMERICAL RESULTS AND DISCUSSIONS 



We will not show the details of Solving the full Salpeter equation, only give the final 
results, interested readers can find the the detail technique in Ref. 

When solving the full Salpeter equation, we choose a phenomenological Cornell potential, 
there are some parameters in this potential including the constitute quark mass and one 
loop running coupling constant. Since we lack data of ++ cc states to determine these 
parameters, we choose the same values as in Ref. [19j which were obtained by fitting the 
mass spectra for 0~ states, but only vary parameter V to fit the mass of XcCb so the following 



4 




FIG. 1: Two-photon annihilation diagrams of the quarkonium. 

parameters are adopted: a = e = 2.7183, a = 0.06 GeV, V = -0.566 GeV, A = 0.2 GeV 2 , 
Aqcd = 0.26 GeV and m c = 1.7553 GeV. With this parameter set, we solve the full Salpeter 
equation and obtain the mass spectra shown in Table I. With the obtained wave function 
and Eq. (j3J), we calculate the two-photon decay width of cc ++ states, the results are shown 
in Table I. To give the numerical analysis of two gluons decay, we need to fix the value of the 
renormalization scale /i in a s (fi). In the case of charmonium, we choose the charm quark 
mass m c as the energy scale and obtain the coupling constant a s {m c ) = 0.36 [ljj. The 
corresponding total decay width are also listed in Table I. 

For the case of bb system, to determine the parameters, we fit the available masses 
M Xb0 = 9.8599 GeV, M x > w = 10.2321 GeV, and we also fit the mass of rj b at 9.364 GeV Q. 
So the parameters are Vq = —0.553 GeV, Aqcd = 0.20 and rrib = 5.13 GeV, other parameters 
are same as in the case of cc. With this set of parameters, the coupling constant at scale 
of bottom quark mass is a s (rrib) = 0.232. The corresponding mass spectra, two-photon and 
total decay widths are shown in Table II. 

TABLE I: Two-photon decay width and total width of P-wave ++ charmonium states, where the 
total width is estimated by the two-gluon decay width Ttot — ^2g- 

Mass GeV T 2l keV Ttot MeV 



Xco(lP) 3.4159 3.78 10.3 

Xco( 2P ) 3 - 8311 3 - 51 9 - 61 

Xco( 3P ) 4.1324 2.86 7.83 

Xco(4P) 4.3694 2.42 6.62 



We compare our predictions with recent other theoretical relativistic calculations and 
experimental results in Table III. Our results of T^ agree with the estimates of Refs. 
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TABLE II: Two-photon decay width and total width of P-wave ++ bottomonium states, where 
the total width is estimated by the two-gluon decay width Ttot — ^2g- 



Mass GeV 



Xbo(lP) 
Xw(4P) 



r 27 eV 



9.8601 
10.2239 
10.4970 
10.7192 



48.8 
50.3 
44.7 
41.9 



tot 



MeV 



0.887 
0.914 
0.813 
0.761 



TABLE III: Recent theoretical and experimental results of two-photon decay width and total width. 

keV r*' t MeV keV eV r** t MeV eV 

Ours 3.78 10.3 3.51 48.8 0.887 50.3 

Crater[3] 3.96, 3.34 

Gupta [8] 6.38 13.44 80 2.15 

Huang [9] 3.72±1.11 12.5±3.2 

Ebert [10] 2.9 1.9 38 29 

Miinz [11] 1.39±0.16 1.11±0.13 24±3 26±2 

E835 [20] 9.8±1.0±0.1 

CLEO [21] 3.76±0.65±0.41±1.69 

PDG [22] 2.87±0.54 10.4±0.7 

BES [23] 12.6l^t?;? 
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211 . All the values of r* c ° listed in the table consist with each other. But there are 
large discrepancy in the results of x c qi Xbo and x'bo from different methods, more theoretical 
calculations and experimental measurements are required for these charnels. 

To show the important of the relativistic corrections, we give a estimate of non-relativistic 
results using the same method and with same parameters. The non-relativistic wave function 
can be written as ¥>i(<7)(l+7o) jf, with this wave function we solve the positive part of Salpeter 
equation, and delete the contribution of negative part. The non-relativistic normalization 
condition is / dq8(fi(q) 2 = 2M(2rr) 3 . When calculate the transition amplitude, we only 
keep the lowest order contribution and delete all the higher order contributions by counting 
the exponential number of \q\, the corresponding results are show in Table IV. One can see 
that the relativistic corrections are important especially for the cc system. 
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TABLE IV: Mass spectra and Two-photon decay width with and without relativistic corrections, 
where the "non-rel" means the non-relativistic results. The decay width is in unit KeV for cc 
system and eV for bb system. 





Mass GeV 


non-rel 




non-rel 


Xco(lP) 


3.4159 


3.3931 


3.78 


5.85 




3.8311 


3.8267 


3.51 


5.47 




4.1324 


4.1543 


2.86 


4.61 


Xd)(4P) 


4.3694 


4.4205 


2.42 


3.94 


XtoO-P) 


9.8601 


9.8471 


48.8 


58.3 


X m {2P) 


10.2239 


10.2129 


50.3 


59.9 


Xw(3P) 


10.4970 


10.4906 


44.7 


54.0 


X^P) 


10.7192 


10.7196 


41.9 


48.3 



We comment that in this work we did not include the QCD radiative correction because 
we focus mainly on the relativistic corrections, though it is no doubt that the QCD correction 
is very important and interesting topic. 

In summary, by solving the relativistic full Salpeter equation with a well denned 
form of wave function, we estimate two-photon decay rates: T(x c o —* 27) = 3.78 keV, 
T{ X ' c0 -> 2 7 ) = 3.51 keV, r{ Xb0 -> 2 7 ) = 48.8 eV and T{x' b0 -> 27) = 50.3 eV, and the total 
decay widths: V tot ( X co) = 10.3 MeV, T tot { X ' c0 ) = 9.61 MeV, T tot { Xb0 ) = 0.887 MeV and 
IWxw) = 0.914 MeV. 
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